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Abstract 

Higgs bosons from an extended Higgs sector, such as a two-Higgs-doublet model, can 
have greatly enhanced coupling to the bottom quark. Producing such a Higgs boson in 
association with a single high-pT- bottom quark via gb hb allows for the suppression 
of backgrounds. Previous studies have instead used gg, qq bbh as the production 
mechanism, which is valid only if both b quarks are at high px- We calculate gb — > hb 
at next-to-leading order in QCD, and find that it is an order of magnitude larger than 
gg, qq bbh at the Fermilab Tevatron and the CERN Large Hadron Collider. This 
production mechanism improves the prospects for the discovery of a Higgs boson with 
enhanced coupling to the b quark. 



1 Introduction 



The Higgs boson couples to fermions with strength nif/v, where v = {■\/2Gf)~^^'^ ~ 246 GeV 
is the vacuum expectation value of the Higgs field. Its Yukawa coupling to bottom quarks 
{mi, ^ 5 GeV) is thus very weak, leading to very small cross sections for associated production 
of the Higgs boson and bottom quarks at the Fermilab Tevatron [Ij and the CERN Large 
Hadron Collider (LHC) [2]. However, this Yukawa coupling could be considerably enhanced 
in extensions of the standard model with more than one Higgs doublet, thereby increasing 
this production cross section j2]. For example, in a two-Higgs-doublet model, the Yukawa 
coupling of some or all of the Higgs bosons (/i", A^, H^) to the bottom quark could be 
enhanced for large values of tan/3 = ^2/^1, where Vi is the vacuum expectation value of the 
Higgs doublet that couples to the bottom quark. 

The dominant subprocess for the production of a Higgs boson via its coupling to bottom 
quarks is 66 — *• /t (Fig. [T]),^ where the 6 quarks reside in the proton sea [21 E]- The 6-quark 
sea is generated from gluons splitting into nearly collinear 66 pairs. When one member of 
the pair initiates a hard-scattering subprocess, its partner tends to remain at low pt and to 
become part of the beam remnant. Hence the final state typically has no high-p^ bottom 
quarks. This subprocess may be useful to discover a Higgs boson for large tan/? in the 
decay mode h t^t~ at the Tevatron and the LHC |H jSj, and h at the LHC 

[SI El El • The decay mode /i — 66 is not distinguishable from the overwhelming background 
gg,qq^ hi. 



Figure 1: Production of the Higgs boson via 66 
quarks in the final state. 



h. There are typically no high-p^- bottom 



If one instead demands that at least one 6 quark be observed at high p^, then the leading- 
order subprocess for associated production of the Higgs boson and bottom quarks is gh hh 
(Fig. 12)) 13 E]-^ The presence of a high-pr bottom quark in the final state has distinct 
phenomenological advantages since it can be tagged with reasonably high efficiency. In the 
case of /i — >• r+r~, /i+yU~ the 6 quark can be used to reduce backgrounds and to identify the 
Higgs-boson production mechanism Cm EI- The trade-off is that the cross section for 
gh hh., with the 6 quark at high p^, is less than that of 66 h. 

If the Higgs boson decays via h 66, the presence of an additional high-p^ bottom 
quark in the final state is essential in order to separate the signal from backgrounds P^ll3j . 
Recent analyses are based on the subprocess gg, qq — > bbh (Fig. and demand a final state 
with four jets, with either at least three 6 tags, or with four 6 tags HH CH CHI UHl HTj . 

^We use h to denote a generic Higgs boson. In a two-Higgs-doublet model, h may denote any of the 
neutral Higgs bosons {hP , , A^). 

^This includes the charge-conjugate subprocess gb hb. All charge-conjugate subprocesses are under- 
stood throughout this paper. 
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Figure 2: Associated production of the Higgs boson and a single high-pj- bottom quark. 



However, the cross section for this subprocess is less than that of gb — > hb. We therefore 
suggest that it may be advantageous to search for h bb hj demanding just three jets in 
the final state, all of which are b tagged [T^ The three-jet final state will have bigger 
backgrounds than the four-jet final state, but the significance of the signal {S/\fB) is likely 
to increase. 

It is only valid to use gg, qq bbh as the production subprocess when both b quarks are 
at high pt- If only one of the b quarks is at high IH E], the integration over the 
momentum of the other b quark yields a factor ln{mh/mf,) which invalidates perturbation 
theory. Our calculation of gb hb sums these logarithms to all orders, and results in a 
well-behaved perturbation series. 




Figure 3: Representative diagrams for associated production of the Higgs boson and two 
high-pT bottom quarks: (a) gg — > bbh (8 diagrams); (b) qq bbh (2 diagrams). 

In this paper we calculate the cross section for the associated production of the Higgs 
boson and a single bottom quark {gb —>■ hb) at next-to-leading order. We provide results 
for both the Tevatron and the LHC The cross section for the subprocess 66 — * /i is already 
known at next-to-leading order The cross section for the subprocess gg,qq — > bbh, 

which has two high-p^ bottom quarks, is known only at leading order, but the analogous 
subprocess gg, qq tih has been calculated at next-to- leading order [1111201, so the next-to- 
leading-order result for gg,qq bbh could be made available. Thus our calculation completes 
the set of next-to-leading-order cross sections for the subprocesses bb h, gb hb, and 
gg, qq bbh. 

In Section |21 we discuss the leading-order cross section for gb — > hb. In Section 01 we 
discuss the correction of order 1/ ln(m/i/mft), due to initial gluons splitting into hb pairs. In 
SectionlHwe discuss the correction of order as] the virtual and real corrections are discussed 
separately. We present our numerical results in Sectional Conclusions are drawn in Section IHl 
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Several appendices follow, in which the analytic results and some of the technical details are 
presented. 

2 Leading order 

The leading-order subprocess for Higgs-boson production in association with a single high-p-p 
bottom quark is shown in Fig.|21 Since the scale of the hard scattering is large compared with 
the 6-quark mass, the b quark is regarded as part of the proton sea j2Il 122 123 1211 • However, 
unlike the light-quark sea, the 6-quark sea is perturbatively calculable. This changes the 
way that one counts powers [31 ES] • If the scale of the hard scattering is fi, the b distribution 
function b{x,fi) is intrinsically of order ^^(/i) ln(yu/mb), in contrast with the light partons, 
which are of order unity. This captures the behavior of the b distribution function at low and 
high values of fi, and interpolates between them. As fi approaches from above, ln(/i/mb) 
vanishes; this reflects the initial condition on the b distribution function, b{x, rrib) = 0. As n 
becomes asymptotically large, ^^(/i) ln(/i/mb) approaches order unity,^ and the b distribution 
function becomes of the same order as the light partons. 

With this counting, the leading-order subprocess gb —>■ hb is of order lia{mh/mb) (times 
the Yukawa coupling), where we have chosen the Higgs-boson mass as the relevant scale. The 
leading-order amplitude may be decomposed into a linear combination of two gauge- invariant 
subamplitudes, 

x = ^!;+^^ (1) 

These subamplitudes are gauge invariant in the sense that they each satisfy the Ward iden- 
tity psfj.A'X ^ = 0, where is the gluon four-momentum. They are related to the two 
independent helicity amplitudes for this subprocess. The explicit form of the subamplitudes 
and the helicity amplitudes are given in Appendix A. 

The spin- and color-averaged cross section for gb hb is 

do-gb^hb 1 Q;g(/z) / ybjjJ,) V mj + u'^ 

dt ~ 24 [ ^ ) St ' ^ ' 

where s, t, u are the usual Mandelstam variables (the first and second diagrams in Fig. |21 
have poles in the t and s channels, respectively), 0:5 (/i) is the MS strong coupling, and yb{lj) 
is the MS Yukawa coupling {yb{^)/\/2 = fnb{^)/v in the standard model, where rfib{^i) is 
the MS mass, and V = (v^Gf)^^/^ ^ 246 GeV). We choose the scale fj, = rrih as our central 
value. It is important to use mb{mh) rather than the pole mass when evaluating the Yukawa 
coupling, as the latter is significantly greater than the former, and would yield an inflated 
cross section."^ The cross section for the charge-conjugate subprocess gb — > Kb is identical. 
The cross section is also identical for the production of a pseudoscalar Higgs boson (A°). 

We neglect the 6-quark mass in Eq. Q and throughout, except in the evaluation of the 
Yukawa coupling. This corresponds to the simplified ACOT scheme |231 1^ . The fe-quark 
mass may be neglected, with no loss of accuracy, in any diagram in which the b quark is an 

■^This can be seen by recalling as{y) ~ 2-^ / / Kqc d)) ■ 

''The evaluation of mf,(m/i) is detailed in Ref. We use rn\j(jni,) = 4.2 GeV as the initial condition j26j . 
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initial-state parton. Terms proportional to the 6-quark mass enter only in the 1/ InijUh/rrib) 
correction. This is discussed at the end of the next section. 
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3 1/ hiimh/rnh) correction 



Consider the subprocess gg, qq bbh, shown in Fig. 121 It is of order a| (times the Yukawa 
couphng). Since the leading-order subprocess gb — > hb is of order a'^\ia{mh/mb), this subpro- 
cess is suppressed by 1/ ln(m/i/m;,) relative to the leading-order subprocess (for rrih ^ mt) 

[3123. 

The helicity amplitudes for this subprocess are given in Appendix B. Integration over 
the phase space of the final-state particles is divergent when the b is coUinear with an 
initial gluon,^ since we use mi, = 0. This coUinear divergence is regulated using modern 
dimensional reduction (DR) j2H] , and absorbed into the b distribution function using a dipole- 
subtraction method as formulated in Ref. ^^.^ This subtraction, together with Gribov- 
Lipatov-Altarelli-Parisi evolution of the parton distribution functions, sums terms of order 
agln^{mh/mb), to all orders in perturbation theory, into the b distribution function [211 
1221 12SI 1211 • This yields a well-behaved perturbation expansion in terms of the parameters 
l/ln(m/i/mfc) and as (the latter to be discuss in Section HI). Our final result is in the MS 
factorization scheme. 

Some fraction of the events from this subprocess yield a final state with two b quarks 
at high pt- In that case the contribution of this subprocess to the total cross section is 
enhanced, since either b can be tagged. If the 6-tagging efficiency is e^,, the probability of 
tagging one or more b quarks when both are at high pt is 2eb{l — eb) + e^. This results in an 
enhancement factor of 2 — e?, relative to subprocesses in which only one b quark is at high px- 
If the Higgs boson decays to bb, the enhancement factor remains 2 — e^, if we demand three 
or more b tags and also demand that two of these tags come from the Higgs-boson decay 
products (so that two 6-tagged jets reconstruct the Higgs-boson mass). 

Since we neglect the 6-quark mass throughout the calculation, we are making an ap- 
proximation. To include the 6-quark mass, one would calculate the diagrams of Fig. [HI with 
a finite quark mass [2S1 1211 12Z|- This would introduce terms of order m\/m\ and ml/p^. 
Hence the only approximation we are making by neglecting the 6-quark mass throughout 
the calculation is of order 1/ \n{mh/mb) x m1/m\ and 1/ ln(m/j/m;,) x ml/p^. 

4 as correction 

In this section we discuss the genuine correction of order as- We divide it into two classes: 
virtual and real. CoUinear divergences are isolated and absorbed into the parton distribution 
functions. Soft divergences cancel between the virtual and real corrections. Both types 
of divergences are regulated using modern dimensional reduction (DR) and are cancelled 
using a dipole-subtraction method, as in the previous section. Our final result is in the MS 
factorization scheme. The 6-quark mass is neglected throughout this section; this introduces 
no approximation [23 1211 12Z] • 

'^This pertains only to gg bbh, which makes a much larger contribution than qq bbh. 
^See Ref. |3J for details on the implementation of this method. 
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Figure 4: Virtual correction to gb hb. External-leg wavefunction renormalization diagrams 
(not shown) vanish in modern dimensional reduction for massless particles. 
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4.1 Virtual correction 



The one-loop correction to the subprocess gb hb is shown in Fig. HJ We calculate in 
d = 4 — 2e dimensions, using modern dimensional reduction (DR). We also used conven- 
tional dimensional regularization (CDR) as a check on our calculation [32], as discussed in 
Appendix E. In DR, the result for the one-loop amplitude is 

^ (-2C(.) + D{s, u) - \d{s, t) + C\u)) 

+Cf {C{s) - D{s, u) - C'{t) - C'{u)) + {s<-^t) 
+ A^^^{Ca-C,), (3) 

where the scalar loop integrals C, C, D are defined in Appendix C, and Cp = (-^c ~l)/2A'"c = 
4/3, Ca = Nc = 3. The one-loop amplitude is proportional to the tree amplitude, A'^, except 
for the last term, which is proportional to one of the two gauge-invariant tree subamplitudes, 
A'^ [Eq. ©], times a finite constant. We checked that this amplitude has the structure of 
infrared (soft and collinear) divergences expected from the dipole-subtraction method (see 
Appendix E). 

The above expression contains ultraviolet divergences. These are cancelled by the renor- 
malization of the strong and Yukawa couplings, as discussed in Appendix E. The ultraviolet 
divergences are also regulated using modern dimensional reduction (DR). The renormaliza- 
tion of the Yukawa coupling with this regulator in the MS renormalization scheme is derived 
in Appendix D. 



At 



A^^ 



4.2 Real correction 

The real correction of 0{as) has several contributions. Fig. Efa) shows the contribution 
from real gluon emission, gb gbh; (b) shows the subprocesses qb — > qbh and qb qbh; 
(c) shows the subprocess bb — > bbh; and (d) shows the subprocess bb bbh. Another real 
correction, gg, qq — > bbh, shown in Fig. El is of 0{1/ \ia{mh/mb)); it is discussed in Section|Hl 
The helicity amplitudes for these subprocesses are given in Appendix B. 

The subprocess bb — *• bbh (and bb —>■ bbh) requires some additional consideration. Since 
there are two b quarks in the initial state, this subprocess is of order ag\n'^{mh/mb), which 
is suppressed relative to the leading-order subprocess by q;| ln(m/i/mb). Thus it is not truly 
a correction of order as- Nevertheless, it is a next-to- leading-order correction in powers of 
as and 1/ ln(m/i/mf,), so it is appropriate to include it in our calculation. Furthermore, 
this subprocess yields two b quarks in the final state. Thus, as discussed in Section |21 this 
subprocess is enhanced by a factor 2 — e?, when both b quarks are at high pt- However, this 
contribution is less than one percent of the leading-order cross section, so this point is moot. 

The subprocess bb — > bbh has a contribution from the diagram shown in Fig. Efd) in 
which a gluon splits into a final-state bb pair. Since we neglect the b mass throughout our 
calculation, this subprocess contains a divergence when the b and b are collinear. In reality, 
the &-quark mass regulates this divergence. To approximate this effect, we restrict the bb 
invariant mass to be greater than 2mb- This correctly captures the dominant, logarithmically- 
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Figure 5: Representative diagrams for subprocesses contributing to the real correction to 
gb hb: (a) gb gbh (8 diagrams); (b) q{q)b — > q{q)bh (2 diagrams); (c) bb — > bbh 
(8 diagrams); bb bbh (8 diagrams). 

enhanced term of order Ins /ml. Since this correction is less than one percent of the leading- 
order cross section, this approximation suffices. 



5 Results 

Figures El ^ El show the cross sections for associated production of the Higgs boson and a 
single bottom quark vs. the Higgs-boson mass at the Tevatron and the LHC These cross 
sections pertain to both a scalar and a pseudoscalar Higgs boson. The Yukawa coupling is 
set to its standard-model value. At the Tevatron, the b jet'^ is required to have a minimum 
Pt of 15 GeV and a rapidity of magnitude less than 2, such that it can be tagged by the 
silicon vertex detector; we refer to this as the tagging region. At the LHC the rapidity 
coverage is taken to be \f]{b)\ < 2.5. Two plots are given for the LHC, one with a minimum 
Pt of 15 GeV (appropriate for low-luminosity running) and one with 30 GeV (appropriate 
for high- luminosity running). Each figure has three curves. The curve labeled crLo(l^) is 
the leading-order cross section, calculated with LO parton distribution functions (CTEQ5L 
|33p and couplings evolved at LO, with the factorization and renormalization scales set to 
/i = rrih.^ The notation indicates that there is only one b quark at high pt. The curve labeled 
o"NLo(l&) is the next-to-leading-order cross section, calculated with NLO parton distribution 
functions (CTEQ5M1) and couplings evolved at NLO, with /i = rrih. Only the subprocesses 
that yield a single b quark in the tagging region are included. Some of the NLO subprocesses 
yield two b quarks in the tagging region; this cross section is labeled (Tnlo(2&) in the figures. 
This cross section is dominated by the subprocess gg bbh, discussed in Section El The 

^Partons within a cone of AR — 0.7 are clustered into a single b jet. 

^The evolution of asip) uses the value of Aqcd corresponding to the parton distribution functions. 
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NLO cross section with one or more b tags is given by o"NLo(l^)e6 + o"NLo(2&)(2efc(l — e;,) +el), 
where eb is the 6-tagging efficiency. As is evident from the figures, the NLO cross section is 
dominated by the subprocesses with a single b quark in the tagging region. 

The NLO correction ranges from 50 — 60% of the LO cross section at the Tevatron for 
rrih = 100 - 200 GeV. At the LHC, the correction ranges from 20 - 40% for pt > 15 GeV, 
and 25 — 45% for > 30 GeV, for ruh = 120 — 500 GeV.^ Most of the correction comes 
from the 0{as) contribution. The (9(1/ ln(m/i/mfc)) contribution is small, less than 10% of 
the LO cross section. Thus the terms we are neglecting by using = throughout the 
calculation, of order 1/ ln(m/j/m{,) x ml/m\ and l/ln(m/i/mb) x ml/p"^, are very small. 

As discussed in the Introduction, recent analyses for the decay h bb use gg,q(j ^ hbb 
as the Higgs-boson production subprocess, and demand a final state with four jets, with 
either at least three b tags, or with four b tags HU El US UHl E] . The cross section with 
at least three b tags (two of which come from the decay products of the Higgs boson) is 
o"NLo(2&)(2eb(l — e^) + e^); with four b tags it is (TNLo(2&)efe. Both of these are an order of 
magnitude less than the cross section with three or more jets, with three or more b tags, 
given by (JNLo{^b)eb + crNLo(2&)(2ef,(l — e?,) + ef). Thus our motivation for carrying out this 
calculation was well founded. 

Similarly, the existing studies of h ^ t^t^, fi^fi^ with at least one b tag use gg, qq bbh 
as the Higgs-boson production subprocess [3 HUl E]. One should instead use the NLO 
calculation of gb hb, since this is a much larger cross section. 

The NLO calculation of the cross section for associated production of the Higgs boson and 
a single b quark gives a more accurate estimate of the cross section than the LO calculation. 
This is evidenced by the fact that the NLO calculation of the cross section is less sensitive 
to the choice of factorization and renormalization scales than the LO calculation. Typical 
examples are shown in Figs. El and El where we plot the LO and NLO cross section vs. the 
common factorization and renormalization scale /i, for rrih = 120 GeV at the Tevatron and the 
LHC (solid curves) }^ In Tables [H and |21 we give the cross section evaluated at = nih as the 
central value (these are the numbers plotted in Figs. El- EI); with uncertainties corresponding 
to /i = m/i/2 (upper uncertainty) and /i = 2mh (lower uncertainty). The scale dependence 
is significantly reduced when going from LO to NLO. Our NLO cross section can be used to 
normalize any future studies that make use of this production mechanism. 

Also shown in Figs. IHl and El is the factorization-scale dependence of the cross section, 
with the renormalization scale fixed to /z = m/j (dashed curves). The factorization-scale 
dependence decreases at NLO, as expected. At the Tevatron, the factorization-scale depen- 
dence is negligible, even at LO. At the LHC, the factorization-scale dependence is greater 
than the dependence on the common factorization and renormalization scales. This indicates 
that there is compensation between the factorization and renormalization scales when the 
two are varied simultaneously. 

The uncertainty in the choice of factorization and renormalization scales yields some 
uncertainty in the NLO cross section. In addition, there is an uncertainty in the cross 
section of about 10% due to the uncertainty in the Yukawa coupling (jnbirfib) = 4.2 ± 0.2), 

^This is the size of the correction for /i = m^. The correction is less for smaller values of /z. 

factor BR{h bb)£l is implicit in the following cross sections. 
^^In this and the following figure, we take e;, = 1 when combining (Tnlo(I^) and crNLo(26) to obtain the 
NLO cross section. However, the NLO cross section is dominated by (Tnlo(16) for any value of et- 
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pp^hb @ Tevatron 
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iTih [GeV] 

Figure 6: Cross section for the associated production of the Higgs boson and a single b quark 
at the Tevatron. The b quark is within the tagging region of the sihcon vertex detector 
{Pt > 15 GeV, \ri\ < 2). The curve labeled (Tlo(1&) is the leading-order cross section, 
evaluated with LO parton distribution functions (CTEQ5L) and couplings evolved at LO, 
evaluated at /i = rrih. The notation indicates that there is only one b quark at high pt- The 
curve labeled ctnlo(1&) is the next-to- leading-order cross section, evaluated with NLO parton 
distribution functions (CTEQ5M1) and couplings evolved at NLO, evaluated at /i = nih. 
Only the subprocesses that yield a single b quark in the tagging region are included. The 
cross section for NLO subprocesses that yield two b quarks in the tagging region is labeled 
o-NLo(2&). 
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pp^ hb @ LHC, high luminosity 
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Figure 8: Same as Fig. but at the LHC, and with a 6-tagging region of pr > 30 GeV, 
\r]\ < 2.5. 
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pp^hb @ Tevatron 
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Figure 9: Cross section for the associated production of the Higgs boson and a single b 
quark vs. the common factorization and renormahzation scale /i, for ruh = 120 GeV at the 
Tevatron (sohd curves). The ratio of the cross section at scale fi to the cross section at scale 
fi = rrih is plotted vs. the ratio of the scales. The next-to-leading-order (NLO) cross section 
is less sensitive to the scale /i than the leading-order (LO) cross section. Also shown is the 
dependence on the factorization scale alone, with the renormahzation scale fixed at /i = uih 
(dashed curves). 
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pp^hb @ LHC 
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Figure 10: Same as Fig. ^ but at the LHC. 
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pp^hb @ Tevatron 
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Figure 11: Higgs-boson distribution for associated production of the Higgs boson and a 
single b quark, for rrih = 120 GeV at the Tevatron. At leading order (LO) the Higgs-boson px 
is balanced against that of the b quark, while at next-to-leading-order (NLO) it is balanced 
against that of the b quark and an additional parton. 
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Table 1: Cross sections (fb) for the associated production of the Higgs boson and a single 
b quark at the Tevatron. The central value corresponds to the choice of factorization and 
renormalization scale /i = m^; these values are plotted in Fig.lHl The uncertainty corresponds 
to varying the scale from fi = m/j/2 to fi = 2mh- The b quark is within the tagging region 
of the silicon vertex detector {pt > 15 GeV, \ri\ < 2). The column labeled crLo(l^) is the 
leading-order cross section, evaluated with LO parton distribution functions (CTEQ5L) and 
couplings evolved at LO. The notation indicates that there is only one b quark at high 
Pt- The column labeled crNLo(l&) is the next-to-leading-order cross section, evaluated with 
NLO parton distribution functions (CTEQ5M1) and couplings evolved at NLO. Only the 
subprocesses that yield a single b quark in the tagging region are included. The cross section 
for NLO subprocesses that yield two b quarks in the tagging region is labeled crNLo(2&). 
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and of about 4% due to the uncertainty in the strong coupling [211 ■ The uncertainty in the 
gluon distribution function (which also reflects itself in the uncertainty in the b distribution 
function) is the source of another 10% uncertainty in the cross section [34]. 

Recall that it is only valid to use gg,qq bbh as the production subprocess when both b 
quarks are at high px- To demonstrate this, we evaluated the cross section for the production 
of the Higgs boson and one high-pj^ b quark using this subprocess by integrating over the 
momentum of the other b quark. For rrih = 120 GeV, this underestimates the NLO cross 
section by a factor of 4.6 at the Tevatron and 2.7 at the LHC. This factor is even larger for 
heavier Higgs bosons. 

We also studied the kinematics of the Higgs boson at NLO vs. LO. The rapidity distribu- 
tion of the Higgs boson remains almost unchanged. The px distribution of the Higgs boson 
does change at low p^, as shown in Fig. ^2 At LO, the px of the Higgs boson is balanced 
against that of the b quark, so the Higgs-boson px cannot be less than the minimum px of 
the b quark. This restriction is lifted at NLO, since the px of the Higgs boson is balanced 
against that of the b quark and an additional parton. 
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Table 2: Same as Table ^ but at the LHC The left side of the table corresponds to a 
6-tagging region of px > 15 GeV, |?7| < 2.5, appropriate for low-luminosity running. These 
cross sections are plotted in Fig. [7| The right side of the table corresponds to pt > 30 GeV, 
\ri\ < 2.5, appropriate for high-luminosity running. These cross sections are plotted in Fig. |S1 
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6 Conclusions 

Previous studies of the associated production of the Higgs boson and a high-pr bottom quark 
have used gg,qq^ bbh as the production mechanism [Tl 1^ 1^ IKj ITTT j ITT ] IT ^ ITH | ITI I ITHl ITT)1 ITTj . 

which is valid only if both b quarks are at high pt- In this paper we have shown that the 
cross section for gb ^ hb is an order of magnitude larger than that of gg,qq — > bbh. 

This production mechanism improves the prospects for the discovery of a Higgs boson with 
enhanced coupling to the b quark. We evaluated the cross section for this subprocess at the 
Tevatron and the LHC at next-to- leading order in QCD. These cross sections can be used to 
normalize any future studies of this production mechanism. They pertain to both a scalar 
and a pseudoscalar Higgs boson. We have included gb — > hb in the multi-purpose NLO 
Monte Carlo program MCFM ^Tj I35j. We encourage studies of the signal and backgrounds 
for associated production of the Higgs boson with a single high pt bottom quark. 
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Appendix A 



In this Appendix we present explicit results for the leading-order subprocess gb hb. In 
order to be systematic, we give results for the (unphysical) crossed subprocess bbgh, in 
which all particles are taken to be outgoing, as shown in Fig. El The amplitudes for the 
physical subprocesses gb — >■ hb and gb — >■ hb may then be obtained by crossing. The 6-quark 
mass is neglected throughout. All expressions are presented in d = 4 — 2e dimensions, using 
modern dimensional reduction. 



b,pi ^ , b,p2 




9,P3 h,P4 



Figure 12: Four-momenta (all outgoing) of the particles for the (unphysical) subprocess 
— >■ bbgh. The arrows indicate the flow of fermion number. 

The amplitude for the leading-order subprocess bbgh may be written in terms of 
the four-momenta of the b, b, and gluon. It is a linear combination of two gauge-invariant 
subamplitudes, 

Ao ^ ej3*(p3)^[; = ej3*(p3)(^!^ + ^^) (4) 
= ^^2e^5l^v^T"2(2'^^|l-'^^) f^-^) (5) 

V2 V'S23 Sis J 

A^B = ^^^''gs^V2T'^^^^^'^{2'W3\l''') , (6) 

V2 -523513 

where yb is the MS Yukawa coupling {yi,{n)/V^ = rni,{jj)/v in the standard model, where 
m(/i) is the MS mass and v = {\/2Gf)^^^'^ ~ 246 GeV), T° are the fundamental-representation 
matrices of 5f/(3) (Tr T"T'' = 5"V2), s^j = (pi+Pj^, Ps = 7/.P3, |1^) = v^{pi), (2±| = 
u^{p2), is the 't Hooft mass (introduced such that the renormalized couplings are dimen- 
sionless in d dimensions), /12 denotes the helicity of the b quark (the b has the same helicity 
in the massless approximation), and h^ denotes the helicity of the gluon. The subamplitudes 
are gauge invariant in the sense that they each satisfy the Ward identity Psf^A'X s = 0. 

One may also describe this subprocess in terms of helicity amplitudes. We define helicity 
amplitudes. A, with the overall factors of the coupling constants removed, 

A(i^,2^,3;=') = /x2^^?5^v^TM(l^^2^^3J^), (?) 

where the Higgs-boson four-momentum is tacit. There are two independent helicity ampli- 
tudes, 



^(1-^2+3 



b^H^^g) '(13) (23) 
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. [12]^ 



[13] [23] 



(9) 



using the spinor inner-product notation as reviewed in Ref. ^36j. The hehcity- reversed am- 
phtudes are equal to these (times —1 if the Higgs boson is a pseudoscalar) , by parity. The 
hehcity amphtudes are related to the amplitudes A'X and A'^ above by 



A(lf,2^3±) 



Squaring the amplitude and summing over colors and helicities gives 



(10) 
(11) 



E \M' = sgn(.,2)16/^^7| ( ^ 



col.hel 



S13S23 



The spin- and color-averaged cross section for the physical subprocess gb 
then obtained from the above by crossing (si2 u, S13 s, S23 t). 
In the dipole-subtraction method, the tensor 



(12) 

hb, Eq. 0, is 



col,h2=±l 



,Sl3 + S23) 
"Sl3"523 



(pi 


_pr\ 


(Pa. 




U23 


S13J 


\S23 


S13J 



mls23 



S13S23 
is also needed. 



S13S23 



»r + P3PI) - r'''' iP3P2 + P3P'2) 

S13S23 



(13) 



Appendix B 

We present the helicity amplitudes for the 2^3 subprocesses shown in Figs. El and El All 
amplitudes are calculated in d = 4 dimensions. The amplitudes in c? = 4 — 2e dimensions 
using modern dimensional reduction may be obtained via gs gsf^'', Vb — ^ The 
calculations were checked with the code MADGRAPH P7j . 

The subprocesses shown in Figs. Efa) and El^a) may be obtained from the helicity ampli- 
tudes for the unphysical subprocess bbggh. These helicity amplitudes may be written 
as 



^(1-^,2^, 353, 4j.) = {gsV2f^^ + 



(14) 



where Ag and Aa are the symmetric and antisymmetric combination of color-ordered ampli- 
tudes. The three independent helicity configurations are given by {sijk = {Pi + Pj + Pk)"^) 



yl,(l5,2,,3,,4J- ^^13^ ^14^ ^23) (24) ^^^^ 



A^(l^,2+;3+4+)=i—- ^( \ (16) 

^ ^' (34) V (13) (24) (14) (23); ^ ' 
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2+;3;,4;; 



[12f 



[13] [14] [23] [24] 



[12]! 1 I 1 ' 
'[34] ^,[13] [24] + [14] [23], 

[12] Si23 



(17) 
(18) 



+ 



V(13) (23) [14] [24] 
[13] ((14) [12] -(34) [23]) , ((24) [12] + (34) [13]) [23] " 



(13) [14] Si34 



(23) [24] s 



(19) 



234 



^a{^l , 2f|"; 3^, Ag 



(34) [23] - (14) [12] / [23] , (-sis + su) [13]' 



(13)S34 V [24] 

(24) [12] + (34) [13] / [13] , (s23 - S24) [23] ^ 



[14] S134 



(23) S34 



.[14] 



[24] S234 



(20) 



Since these are not color-ordered amplitudes, the order of the gluons in the arguments of the 
functions Ag and Aa is irrelevant. Squaring and summing over colors gives 



El^(l-^,2^,3J^4J^)p 



col 



4 / Vb 



-I N? 



Nr. 



lAJ'^ + NAAal 



(21) 



The four-quark subprocess of Figs. Efb) and Efb) may be obtained from the helicity 
amplitudes for the unphysical subprocess — * bbqqh, 



A (-I hi nh2 nhs a hi 



Vb 1 



The two independent helicity configurations are given by 

A(l-+ 2+ 3+ 4,-) = /(I, 2, 3, 4) + /(2, 1, 3, 4) 
A(l-+ 2+ 3r, 4+) = /(I, 2, 4, 3) + /(2, 1, 4, 3) 



where 



/(l,2,3,4) = z 



J13] ((14) [12] -(34) [23]) 



Squaring and summing over colors gives 



■S34'5l34 



^ l^(l-^2^,3J^4J^)p = g's (^^j {N^ - 1)1^ 



(22) 

(23) 
(24) 

(25) 
(26) 



The four-6-quark subprocesses in Figs. Efc) and (d) may be obtained from the helicity 
amplitudes for the unphysical subprocess bbbbh. These helicity amplitudes may be 
obtained from the four-quark amplitudes above by subtracting a term with one pair of the 
identical quarks exchanged, 



^(l-^^2^^3-^^4^^) 



Nr 



■'12 «4 



Nr 



(27) 
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where 



A = A(l^^2;:^3-^^4^) (28) 

The four independent hehcity configurations are given by 



M^t^ 2fJ", 3^ 


,4,-) 


= /(l,2,3,4) + /(2,l,3,4) 


(30) 


^(Ifc 5 2j, , 3^ 




= /(l,2,4,3) + /(2, 1,4,3) 


(31) 




,4^) 


= /(3,4,l,2) + /(4,3,l,2) 


(32) 


^(Ife 5 2^ , 3^ 


,4^) 


= /(3,4,2,l) + /(4,3,2,l). 


(33) 



Squaring and summing over colors gives 

^I^(1-^2^,3J^4J^)P = ^;^(^^^ iN!-l)(\A\' + \A^^\' + ^Re{A^^A*)^ . (34) 

Appendix C 

In this appendix we hst the scalar integrals which result from the Passarino-Veltman tensor 
reduction of the one- loop diagrams of Fig. HI as given in Eq. Q. We define {d = 4 — 2e) 

ri I 2 2 2 2 2 2\ 

^ J (27r)^ [P - ml][{k + pi)2 - ml][{k + p, + P2Y - ml] 

n/'22222 2 2 2 2 2\_ 
-Do (Pl , P2 > 5 P4 , Pl2 > P23 ; "^0 > "^1 5 "^2 5 "^3 ) = 

2e [ ^'^k 1 

^ J {2txY [P - ml] \{k + piY- mf] [(fc +pi+ pa)^ - m^] [(A; + P1+P2+ p-if - ml] 



where pf, = {pi + PjY, and 



The scalar integrals needed are [s and t are generic invariants here) 

5o(0;0,0) = (39) 
^0(^^0,0) = ^(; + 2-ln^) .<0 (40) 

Co(0,0,.; 0,0,0)^^ = ^2-(-2--^^^ + W^] ^<0 (41) 
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Co(m^,0,s; 0,0,0) = 



Do(0,0,0,m2,s,t;0,0,0,0) = 



mi-s 



IGtt^ m\ — s 
s,ml<0 



iD{s,t) 

7t 
icr ( 



- - In" 



— S 



-In 



e -mt 



2 



167r2 
1 



-mt 



+ 



In^ 



— mt 



2 

~t 

+ ln' 



1 1 

t 



In 77 + m ^ 



-mt 



-mt 



R 



— m 



ft/ 



2 ' 2 



(42) 



(43) 



m^j < s, t < 



R(a;,|/) = Inxlny — lnxln(l — x) — ln?/ln(l ~ ?/) + — — ^h{y) 



(44) 



The analytic continuation of the above results to the physical region is accomplished through 
the use of 

ln(— s — irj) = In \s\ — inQ{s) , (45) 

where s is a generic invariant (including m^) and 1] is a small positive number and, for the 
dilogarithms with arguments greater than unity, by means of 



Re [Li2(x)] = -Lia - H In^ a; . 



X 



TT 



(46) 



Appendix D 

In this appendix we derive the QCD counterterm for the Yukawa coupling in the MS renor- 
malization scheme jSH]. We consider two different regularization schemes (RS): conventional 
dimensional regularization (CDR) and modern dimensional reduction (DR). 

The Yukawa coupling and the quark mass arise from a common term in the Lagrangian, 

£ = -^QQ{h + v), (47) 

where ?/o is the bare Yukawa coupling, h is the physical Higgs-boson field, and v = (v^G^?)"^/^ 
~ 246 GeV is the vacuum-expectation value of the Higgs-doublet field. It is evident from 
Eq. ()47p that the bare quark mass is related to the bare Yukawa coupling by mo = yov / V2. 

The quark mass and Yukawa coupling receive corrections at one loop in QCD. We express 
the bare parameters in terms of the MS values and a counterterm, 

yo = fi'y{l + Syns) (48) 
mo=m + 5mRs, (49) 
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where the subscript on the counterterm indicates that it depends on the regularization scheme 
(RS). The parameter /i is the 't Hooft mass, introduced to keep the renormahzed Yukawa 
couphng dimensionless in d = 4 — 2e dimensions. The Higgs vacuum-expectation value does 
not receive a correction at one loop in QCD. From the above equations and the relation 
^0 = Z/o'^^/v^ we find that the mass and Yukawa-coupling counterterms are related by 
^Urs = ^^Rs/^- Thus we may obtain the Yukawa-coupling counterterm from the mass 
counterterm. 

The one-loop quark propagator is given by 

(50) 



p-mo + Srs(p) p-m- (5mRs + T.rs{p) ' 

where zSrs(p) is the one-loop quark self energy. Since it is ultraviolet divergent, it depends 
on the regularization scheme. The position of the pole in the propagator at one loop is^^ 

mpoie = m + 5mRs - SRs(m) . (51) 

The MS mass is defined via 

5"^CDR = S(m)|div , (52) 

where S(m)|div is the divergent part, proportional to cr/e, of the quark self energy (which 
is the same in CDR and DR). Since the pole mass is a physical quantity,^^ independent of 
the regularization scheme, Eqs. and yield 

SmY)K = S(m)|div + SDR(m) - Scdr("^) • (53) 
The one-loop quark self energy (in 't Hooft-Feynman gauge) is 



zSrs(p) = [ -—^{zgsfi'Yn—^ (WVT" 

J (Zirr r)-\-k-'m. 



( [(-2 + 2e)p + (4 - 2e)m] B{p'^) + (-2 + 2e)pA{p'^) in CDR 

= -glcA (54) 

[ (-2p + 4m)S(p2) - 2pA(p2) inDR 

where Cp = (N^ - l)/2N^ = 4/3 and 



{2tiY P [{p + kf - 



m 



B{p^) = Bo{p^;0,m^) . (56) 

Using 



^^At one loop, the mass m in the argument of the quark self energy may be regarded as the pole mass or 
the MS mass. 

^•^The quark pole mass is a physical quantity within perturbation theory, which suffices for our purposes. 
However, it is unphysical once nonperturbative QCD is taken into account |39[ I40| . 
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gives 



Air 





[3 , , 1 




- + 4 + 5rs 



m 



(59) 



SRs("^) ■ 

where 5cdr = and (5dr = 1- The counterterm is then obtained from Eqs. and 

Sturs = -^Cpcr -+Srs rn . (60) 
Using 5?/Rs = 5"^Rs/^ then gives 



RS 



47r 



-+5 



RS 



(61) 



This resuh is independent of the gauge chosen. 

The MS counterterm for the strong couphng, g^g = fi''gs{l + Sgg-^), analogous to Eq. (jHU, 

is 



x RS '^•5 



h — dRs 

e 6 



(62) 



3, Tp = 1/2, is the number of 



where bo = Po/2 = (11/6)C^ - (2/3)Tpn/, Ca 
hght quarks, and 5cdr = 0, 5dr = 1- 

The relation between the pole mass and MS mass may also be obtained from Eqs. ()51|). 
(jS21), and (jSHD, 



m 



pole 



Air 



(63) 



Appendix E 



In this appendix we discuss two of the checks performed on our calculation of the virtual 
correction presented in Section 14.11 and Appendix C. We checked that the structure of the 
infrared (soft and collinear) divergences is as expected from the dipole-subtraction method 
PU|. We also performed the calculation in both conventional dimensional regularization 
(CDR) and modern dimensional reduction (DR) and verified the scheme independence of 
our results. 

The structure of the divergences for gh hh at one loop is 



2Re {AiAD 



div 



I A la^S" 



m 



X 



--{Ca + 2C, 



h/ 

In A + In -4 



m 



h 



m 



- {Ca - 2Cf) In ■ 



-u 



h, 



(64) 



where Cp = {N^ — l)/2Nc = 4/3, Ca = Nc = 3 and where is the tree amplitude given in 
Eqs. dl)) and ((Zj), and Ai = e|j^*(p3)^i where Ai is the one-loop amplitude given in Eq. (jHl). 
This expression contains both infrared and ultraviolet divergences. The latter are removed 
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by renormalizing the Yukawa and strong couplings using the counterterms given in Eqs. 
and respectively. This leaves the infrared-divergent expression 

2Re{A^A*) + \Ao\'2{Syns + Sgf) 

+ - Ica (in ^ + In ^ ) - {Ca - 2Cf) In ^ - 6o - 

which has the structure expected from the dipole-subtraction method [20] • 

We also calculate the relation between the virtual amplitude in CDR and DR. We find 

2Re (^1^*) IcDR = 2Re (^i^*) |dr - \Ao\^^2Cf . (66) 

The above relation is consistent with the set of rules given in Ref. j2H|, augmented by the 
regularization-scheme-dependent renormalization of the Yukawa coupling given in Eq. 



div 



--(Ca + 2Cf) 
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